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Abstract 

We analyze resonant tunneling subject to a non-adiabatic time-dependent bias-voltage through 
an asymmetric single molecular quantum dot with coupling between the electronic and vibrational 
degrees of freedom using a Tien- Gordon-type rate equation. Our results clearly exhibit the ap- 
pearance of photon-assisted satellites in the current- volt age characteristics and the elimination of 
hot-phonon-induced negative differential conductance with increasing ac driving amplitude for an 
asymmetric system. This can be ascribed to an ac-induced suppression of unequilibrated (hot) 
phonons in an asymmetric system. 

PACS numbers: 85.65.+h, 71.38.-k, 72.40. +w, 73.63.Kv 
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Recently, phonon-mediated dc tunneling through a quantum dot (QD) with cou- 
pling to an internal vibrational (phonon) mode (IVM) has attracted much interest, both 
experiment ali'^'^ and theoretical.-i^>^i^i^i2i>i2iiii Electronic transport measurements in sus- 
pended Carbon nanotubes (CNT) indicate that the current-voltage characteristics display 
perfect signatures of phonon-mediated tunneling, e.g. stepwise structures having equal 
widths in voltage and gradual height reduction by the Franck-Condon (FC) factor.-!^ Fur- 
thermore, striking negative differential conductance (NDC) and current peaks have been 
observed at the onsets of phonon steps,- which are ascribed to the combination of highly 
asymmetric tunnel-coupling rates and the voltage-triggered unequilibrated (hot) phonon 
effect .-iiSiiii On the other hand, ac-transport through a single molecule has also attracted 
considerable interest, but without consideration of electron-phonon coupling (EPC).— In this 
letter, we will examine ac-transport in the presence of EPC when a time-dependent bias- 
voltage is applied between the two electrodes in the non-adiabatic regime.— We find that 
application of a strong ac-amplitude causes suppression of hot phonons and thus elimination 
of NDC for a large asymmetric system. 

We consider a quantum dot (QD) with one spinless level coupled to two electrodes, and 
also linearly coupled to a molecular IVM. The model Hamiltonian is 

H = Hleads + Hjnol + Hb + Ht, (la) 

with 

Hleads = J^^^kcJkC^k' (lb) 

Hmol = £dc\c^ + wofl^a + ^c\c^{a) + a), (Ic) 
Ht = 5^(VikCd + H.c.), (Id) 

r),k 

where cj^^ (^rjk) is the creation (annihilation) operator of an electron with momentum k 
in lead r] {rj = L, R). The energies Enkit) = + ^U.r,{t) include rigid shifts of the Fermi 
energies of the electrons in the leads due to the applied time- dependent bias- voltage. Unit) = 
+ M^cos(fit), with f/° (ujj) being the dc(ac) part of bias- voltage and Q is the driving 
frequency. c);(q) is the creation (annihilation ) operator for a spinless electron in the QD. 
a^{a) is the phonon creation (annihilation) operator for the IVM with energy Uq. A represents 
the EPC constant; describes the tunnel- coupling between the QD and lead rj. Here, we 



assume that the Fermi energies of two leads are zero at equihbrium and [/^ = —U^ = eV/2, 
ul = —ur = eVac (V and Vac are the dc bias- voltage and ac-bias amplitudes, respectively). 
In the following, we will use units where h = ks = e = 1. 

It is well-known that the electron-phonon interaction term in Eq. flTc]) can be eliminated 
by a canonical transformation,— leading to a renormalization of the parameters, Sd = Sd — gX 
{g = \/ujq), and of the tunnel-coupling, V^exp[(yf(a^ + a)]. In the weak tunneling regime and 
high temperature approximation, <ti T, ujq (F^ is the tunneling rate of lead r] and T is 
the temperature), rate equations for the electron-phonon joint probabilities, Pqq and p"^ for 
zero- and one-electron together with n excited phonons on the molecule (incorporated with 
the FC-modified tunneling rates) are physically appropriate for the description of resonant 
tunneling through a single molecule involving EPC-'^'^iSi'i^i'ii On the other hand, in the limit 
of high driving frequency, ^ F, T, of interest in this letter, the ac-bias oscillates so fast 
that an electron experiences many cycles of the ac-bias during its presence inside the dot, 
and thus can not sense the details of the dynamics within a single period. Correspondingly, 
the rate equations can be established by directly applying a Tien-Gordon-type tunneling 
rate in the presence of such an ac-bias in the non-adiabatic limit,— 

nm Poo), (2) 

m 

Pll — ''^^(XmnPoO ~ ^mnPll)y (3) 
m 

with the normalization relation YlniPoo~^ Pii) ~ ^- "^^^ electronic tunneling rates are defined 
as 



oo 

^nm — '^r],nm — '^ri'Jnm ^Jj (^'q) 
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Tj r) ]=-co 

xf^{ed+ {m-n)uJo- jQ), (4) 

oo 2 
rj r) j=-oo 

X [1 - fr,{ed + (m - n)uJo - J^l)], (5) 

with the FC factor given by {p = min{m, n} and q = max{m, n}, denoting the smaller and 
larger of the quantities m and n, respectively) 



7n^ = e-^^/l"-"l^[Ll^-"l(/)]^ (6) 
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where L'^ix) is the generahzed Laguerre polynomial. Jj{x) is the Bessel function of order 
j. frji^) = [1 + e^^~^^^^'^]~^ is the Fermi-distribution function. Moreover, the dc current can 
be calculated by solving the rate equations, Eqs. ([2]) and ([3]), in the steady state condition 
and evaluating the net tunneling rate through one of electrodes (for example, the left lead): 

^ = E(rL™Poo-rz,^nPn)- (7) 

nm 

Note that if the driving amplitude vanishes, = 0, the above equations ©-([Tj) reduce 
exactly to the rate equations formulated without an ac-bias.-'^'^'^i'i^i'ii 

In our numerical calculation, we set the parameters as: cuq = 1 is the energy unit, 
^r/^l = 10^, g = I, Ed = 0, T = 0.02ci;o5 and the driving frequency is = 0.5uq. Note 
that we choose a large asymmetry in the tunneling rates, F/j/Fi = 10^, and an intermediate 
EPC strength, g. The transport properties of this specific system exhibit weak NDC at the 
onsets of phonon steps under hot phonon condition.— In particular, numerical fits of the 
experimentally measured data for the I-V curves show (7 ~ 1 for long CNTs.- 

Figure 1 shows the dc current-bias- voltage characteristics as a function of various driving 
amplitudes. For comparison, we also plot the results without an ac-bias. Vac = 0, as the 
solid line, showing (1) step behavior with width equal to 2uJo, (2) the FC factor-suppressed 
height, and (3) NDC at the onsets of phonon steps starting from the second step. When an 
ac-bias is applied, it is obvious to see from Fig. 1 that a photon- assisted tunneling feature — 
constituted of new steps atV = 2nQ {n is integer) — is superimposed on the phonon-mediated 
I-V characteristics. For stronger driving amplitudes, the suppression of current is more 
pronounced in the weak dc-bias region, V < 2Q. However, in the strong dc-bias region, 
V > Aujq, the application of the ac-bias always enhances the current. Such nonuniform 
behavior is reflected in the structure of the Bessel function in the tunneling rates, Eqs. (jl]) 
and ([5]). More interestingly, we observe that for a weak ac-amplitude. Vac = 0.2^010, the 
position of NDC moves to the ends of the phonon steps, while in the cases of strong ac- 
amplitudes. Vac = O.Scuq and I.Ouq, the NDC even disappears. It should be noted that strong 
vibrational relaxation can result in the elimination of NDC even in the absence of an ac- 
bias.'iS''-^ In that case, photon-assisted steps still remain in the current-voltage characteristics 
when an ac-bias is applied (not shown here).— 



To gain insight into this peculiar behavior, we adopt the method of Ref. [id by examining 



the phonon occupation numbers (PONs), Pqq. For the positive dc bias- voltage range of 



interest in this letter, > 0, the QD is always in an empty occupation state owing to 
Tr ^ Tl (so that the electron has a much stronger tunneling-out rate to the right lead 
than the tunneling-in rate from the left lead). Therefore, a nonzero value of indicates an 
open channel contributing to the current. In this case, we can rewrite the current formula, 
Eq. (j7]), with good accuracy, as 

N 

I = J2Voo, (8) 

n 

with (A^ denotes the maximum number of PONs of the open channels) 

m m j=—oo 

xQ{-Id- {m-n)uJo+jn + V/2) (9) 

[9(x) is the Heaviside step function]. Starting from these two equations, we interpret the 
NDC and positive differential conductance (PDC) behaviors of the I-V characteristics in 
the presence of the ac-bias in Fig. 1. 

Figure 2 displays the Pqq under several dc-bias voltages and ac-amplitudes as well. It is 
evident that for small dc-bias, Vi = 2.0ujq (the panels of the first row in Fig. 2), only the 
channel with the PON n = is active, Pqq ~ 1, and the application of an ac-bias has very 
little influence on the PON. In these situations, we have = and thus 

/l ^r^J^TOmJm, (10) 

m 

with J'm = ^^2(m-i) [^i('^)]^- -'■f ac-bias vanishes, the current becomes (we use / here 
to denote the current in the absence of an ac-bias) 

Ji ~ rL(7oo + 7oi). (11) 

When the dc-bias exceeds the transition point Vi = 2.0co'o, for example V2 = 2.3^0 (the 
panels in the second row of Fig. 2), two channels, n = and 1, are open in the absence of an 
ac-bias, which indicates an unequilibrated (hot) phonon distribution. Correspondingly, the 
current is 

h ^ ri[Poo(7oo + 7oi) +pio(7io + 7ii +712)] 

- ^Lloo + rLPoo(7ii + 7i2 - 7oo), (12) 
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because Pqq + pjg ~ 1. As a result, I2 — h — TlpIqC'^ plQg^{g^/2 — 1) < due to g = 1 
in this letter, leading to the appearance of NDC— When an ac-bias is applied, we find — 
strikingly — from Fig. 2 that the hot phonon is suppressed and the n = 1 channel becomes 
closed. This suppression of pjg stems from a weak increase of p"^ due to photon-assisted 
processes. As a result, we have I2 ~ /i and the NDC reverts to a pure PDC at the onsets of 
the phonon steps. Furthermore, for V3 = S.Ocuq and V4 = S.Suq, two channels, n = and 1, 
are activated by the dc-bias and we readily obtain I3 ~ 14. The ac-bias-induced suppression 
of hot phonons changes the current of the case of V3 as /a ^ Ji. However, the application 
of a weak ac- amplitude. Vac = 0.2070, can not completely suppress the n = 1 hot phonon at 
the moderately stronger dc-bias V4, and thus 

m m 

= h + TLplo^illmJL-lOmJm), (13) 
m 

with = ^°l2m-5['^i('^)]^- "^^^ second term of the last line of Eq. (fT3l) may be positive 
or negative depending on the EPC constant g and the ac-amplitude Vac- For the weak 
ac-amplitude Vac = 0.2u;o, it is negative, leading to NDC; while for stronger amplitudes, 
Vac = 0.5ujq and I.Oujq, it is positive and PDC is observed. 

In summary, we have employed the Tien-Gordon-type rate equations incorporated with 
the FC model to analyze photon-assisted resonant tunneling through an asymmetric single 
molecular QD with EPC in the non-adiabatic and high frequency regime. Our results predict 
that the application of an ac-bias causes suppression of the dc-bias-triggered hot phonon, 
leading to the elimination of NDC in the I-V characteristics. 

This work was supported by Projects of the National Science Foundation of China, the 
Shanghai Municipal Commission of Science and Technology, the Shanghai Pujiang Program, 
and Program for New Century Excellent Talents in University (NCET). 



6 



^ H. Park, J. Park, A. Lim, E. Anderson, A. Alivisatos, and P. McEuen, Nature 407, 57 (2000). 
2 B.J. LeRoy, S.G. Lemay, J. Kong, and C. Dekker, Nature 432, 371 (2004); B.J. LeRoy, J. Kong, 

V.K. Pahilwani, C. Dekker, and S.G. Lemay, Phys. Rev. B 72, 75413 (2005). 
^ S. Sapmaz, P. Jarillo-Herrero, Ya.M. Blanter, C. Dekker, and H.S.J, van der Zant, Phys. Rev. 

Lett. 96, 26801 (2006); S. Sapmaz, P. Jarillo-Herrero, Ya.M. Blanter, and H.S.J, van der Zant, 

New J. Phys. 7, 243 (2005). 

D. Boese and H. Schoeller, Europhys. Lett. 54, 668 (2001). 

5 A.S. Alexandrov and A.M. Bratkovsky, Phys. Rev. B 67, 235312 (2003); K.D. McCarthy, N. 
Prokof'ev, and M.T. Tuominen, Phys. Rev. B 67, 245415 (2003). 

6 A. Mitra, I. Aleiner, and A.J. Millis, Phys. Rev. B 69, 245302 (2004). 

7 J. Koch and F. von Oppen, Phys. Rev. Lett. 94, 206804 (2005); J. Koch, M.E. Raikh, and F. 
von Oppen, Phys. Rev. Lett. 95, 56801 (2005). 



^ K.C. Nowack and M.R. Wegewijs, |cond-mat/0506552| (2005); M.R. Wegewijs, K.C. Nowack, 

New J. Phys. 7, 239 (2005). 
9 J. Koch and F. von Oppen, Phys. Rev. B 72, 113308 (2005). 

10 A. Zazunov, D. Feinberg, and T. Martin, Phys. Rev. B 73, 115405 (2006). 

11 X. Y. Shen, Bing Dong, X.L. Lei, and N.J.M. Horing, unpublished. 

12 S. Camalet, S. Kohler, and P. Hanggi, Phys. Rev. B 70, 155326 (2004); P. Hanggi, S. Kohler, J. 
Lehmann, and M. Strass, Phys. Rep. 406, 379 (2005); P. Hanggi, S. Kohler, J. Lehmann, and 
M. Strass, Lect. Notes Phys. 680, 55 (2005). 

13 Bing Dong, H.L. Cui, and X.L. Lei, Phys. Rev. B 69, 205315 (2004). 

1^ G.D. Mahan, Many-Particle Physics. (Third edition, Kluwer Academic/Plenum Publisher, New 
York, 2000). 

15 P.K. Tien and J.R. Gordon, Phys. Rev. 129, 647 (1963); C. Bruder and H. Schoeller, Phys. 

Rev. Lett. 72, 1076 (1994). 
1^ W.G. van der Wiel, T.H. Oosterkamp, S. De Franceschi, C.J. P.M. Harmans, L.P. Kouwenhoven, 



cond-mat/9904359 



7 



Figure Caption 



FIG.l: I-V curves as functions of ac bias-voltagc amplitudes, Vac/ujQ = 0, 0.2, 0.5, and 1.0, 
for a fixed driving frequency, Q/u!o = 0.5. The parameters we use in the calculation are: 
5=1, Tl/Tr = 10-3, and T/uo = 0.02. 

FIG. 2: Phonon probability distributions for increasing ac-amplitudc, Vac, under different dc 
bias- voltages, V/ujq = 2.0, 2.3, 3.0, and 3.3. Other parameters are the same as in Fig. 1. 
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Figure 1: I-V curves as functions of ac bias-voltage amplitudes, Vac/wo = 0, 0.2, 0.5, and 
for a fixed driving frequency, 0,/ujq = 0.5. The parameters we use in the calculation are: g 
TL/Tn = 10-3, and T/lvq = 0.02. 
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Figure 2: Phonon probability distributions for increasing ac-amplitude, Vac, under different dc 
bias-voltages, V/loq = 2.0, 2.3, 3.0, and 3.3. Other parameters are the same as in Fig. 1. 
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